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ABSTRACT 



Mean humidity, temperature, wind data and fluctuating 
temperature and velocities data obtained from the R/V Acania 

were interpreted on the basis of empirically derived expres- 

. . 2 
sions describing the temperature structure parameter, C^, , 

and the dissipation rate of turbulent kinetic energy, e . 

2 

Properties of C,p and e are related to a stability parameter 

Z/L and the Richardson number, Ri. Mon-dimensional tempera- 

2 4/3 30 2 

ture structure functions parameters, C™ /Z (-f— ) , and tur- 

1 d Z 

-3 

bulent kinetic energy dissipation rates, Ks ZU... , versus Ri 

were examined for open ocean conditions. 

This experiment verified earlier hypotheses that the pro- 
files obtained from empirically derived expressions using 
overwater data are very much different from those obtained 
using overland data. The distribution of the non-dimensional 
turbulent kinetic energy dissipation rate agreed with that 
predicted by similarity theory. 
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I. INTRODUCTION 



The marine boundary layer is that region in the first 
kilometer over the sea which is defined principally by tur- 
bulent transfer of momentum and heat. In the past decade 
there have been concentrated efforts to make measurements 
in this region and to define physical expressions describing 
the turbulent intensities. The marine environment presents 
a number of problems, especially in accessibility, platform 
stability, and sensor protection. Most important of these 
problems are the lack of an extensive data bases and the 
fact that available turbulence expressions, which are empiri- 
cal in nature, have been formulated primarily from overland 
data. 

With the development of systems dependent on laser tech- 
nology, the Department of Defense interest in turbulence in 
this region pertains to assessing and predicting the environ- 
mental influence on optical wave propagation. The refractive 
index problem for the marine boundary layer is more compli- 
cated than the similar problem over land because of the dif- 
ficulties in making point measurements of the fluctuating air 
temperature from which refractive index statistics for optical 
waves are often obtained. For optical waves the atmospheric 
index of refraction variations primary result from tempera- 
ture-induced density changes. These variations cause scintil- 
lation and phase fluctuations of laser beams propagated in the 
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atmosphere and the loss of sight of objects viewed through 
high magnification telescopes. 

An interdisciplinary group at the Naval Postgraduate 
School has been involved in overwater experiments to relate 
the optical properties and the optically relevant meteoro- 
logical parameters of the marine boundary layer to observ- 
able bulk atmospheric parameters. 

The primary purpose of this work was to evaluate pre- 
viously formulated overland expressions for the overwater 

regime. Expressions presented by Wyngaard et al. (1971) and 

2 . . 

Friehe (19 76) were used to relate C^, to a stability para- 
meter, either the Richardson number or Z/L where L is Monin- 
Obukhov Length. In addition to this, expressions described 
by Busch and Panofsky (1968) were used to relate e to the 
stability parameters. 

Ten minute time averaged values of wind, temperature, 
and relative humidity measured at four levels (4.2, 6.6, 

11.3, and 17.7 meters) were examined in this study. The 
most recent and the majority of the observations examined 
were made during the months of September and October 1976 off 
the coast of San Diego, California. Data previously evalua- 
ted by Hughes (1976) and Karach (1976) were incorporated into 
the final evaluations to supplement the data base. The latter 
data were also measured at four levels (4.2, 6.6, 7.6, and 
13 . 9 meters ) . 
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II. THEORETICAL BACKGROUND 



A. GENERAL 
2 

Cm and e are optically relevant turbulence parameters be 
cause they define the intensity of fluctuations in the refrac 
tive index and the size or scale of refractive index inhomo- 
geneities, respectively. The scales of the temperature 
inhomogeneities that distort image resolution and degrade 
laser beam propagation in the atmospheric surface layer lie 
in the inertial subrange where fluctuations are isotropic. 

The latter characteristic allows theories by Tatarski 
(1961) describing electromagnetic propagation to be applied 
in analyses. Tatarski considered the refractive-index-struc- 
ture parameter, 

0 2 _ [n(x) - n(x+r)] 2 , n ^ 

n ^271 



which is related to the temperature structure function para- 
2 

meter, C T , as 

C 2 = (79 x 10 -6 ^y) -2 C 2 (2) 

i m n 



where P is the pressure and T the temperature of the atmos- 
2 

phere. C^, may be evaluated in a number of methods. 

The dissipation rate e defines the micro-scale, Z Q , as 



Z 



o 




(3) 
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where y is the kinematic molecular viscosity. The micro- 
scale appears in optical parameters for image resolution. 
For example, the exp (-2) folding distance of the mutual 

coherence, p is defined as 
3 o 



P 



o 



1 

(1.5C l ~ l/6 K) 
n o 



(4) 



where K is the average wave number. 

B. TEMPERATURE STRUCTURE-FUNCTION PARAMETER, C T 2 

The temperature structure-function parameter may be 
evaluated or estimated by several different methods depend- 
ing on the measurement procedures. At separations r of the 
order of inertial subrange scales, the temperature structure 
function in a locally isotropic field has the form 

^ 2 _ <[T(x) - T(x+r ) ] 2 > 

C T ^271 (5) 



2 

C^ can also be related to the one-dimensional tempera- 
ture spectrum (The Fourier transform of the correlation func 
tion with separations in the streamwise direction) in the 
inertial subrange as 

S t (K x ) = .25 C T 2 K- 573 (6) 




is the streamwise component of wavenumber. 
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A third approach is based on the hypothesis of Corrsin 
(1951) that in the inertial subrange the one dimensional tem- 
perature spectrum can be expressed as 

S t <K) = axs _1/ 3 K~ 5/3 (7) 



where x is the rate of molecular destruction of <(T-T)> /0 

3. V / Z 

r\ -| 

(in 0 C Z S -X ), e is the rate of molecular dissipation of turbu- 

2-3 

lent kinetic energy (in cm s ) , and a is hypothesized to be 

a universal, dimensionless constant; the same in all flows. 

2 

Equations (6) and (7) yield an expression for C^, m terms of 
the rates of dissipation of turbulent kinetic energy and 
temperature variance, viz. 



C 



T 



2 



$xe 



-1/3 



( 8 ) 



where 8 is an empirical constant with a value of 3.20. 

C. MICROSCALE £ (OR e) 

o 

Both s and x in Equation (8) represent molecular processes 
and forces which act only over very small distances. As pre- 
viously indicated, the particular scale length of interest 
in propagation analyses is the microscale of turbulence de- 
fined by Equation (6) and is dependent on only one dynamic 
parameter, the dissipation rate of turbulent kinetic energy, 

£ , which can be estimated from the one dimensional spectral 
density of velocity as 

S U (K) = E s -2/3 K- 573 (9) 

where E is an empirical constant determined to be 0.51. 
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D. TAYLOR'S HYPOTHESIS 

Because measurements are normally made at a point a rela- 
tion between temporal (f) and wavenumber or spatial frequency 
(k) is required. According to Lin (1953), if one considers 
that eddies of wavenumber K are correlated only over distances 
not large compared to 2Tr/k, then over such distances the non- 
uniformity of convection velocity should be small relative 
to the convection velocities. Further, if the mean velocity 
is a good approximation to the convection velocity, then 
0 >> U' where O' is the spatial gradient of U. If these cri- 
teria are met the turbulence field may be considered as being 
frozen and swept at the mean wind speed 0 past a stationary 
instrument. This is known as Taylor's hypothesis and it 

allows frequency spectra to be interpreted as streamwise 

2 

wavenumber spectra. This enables the evaluation of C,p from 
Equation (6) using point temperature fluctuation measurements. 

E. RELATING C T ? AND e TO STABILITY PARAMETERS 

Turbulence similarity theory relates both e and x to 
. . 2 

measurable stability parameters so that C,p can be estimated 
using Equation (8). The 1968 Air Force Cambridge Research 
Laboratory (AFCRL) boundary layer study yielded information 
on the behavior of both e and x and the combination of dimen- 
sional and physical arguments that make up surface-layer 
similarity theory. The AFCRL results demonstrated that the 
statistics of both the mean and turbulent fields, when pro- 
perly non-dimensionalized , are universal functions of the 
stability parameter (Z/L). 
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Non-dimensionalized e and x values can be expressed as 



— , = f -i (Z/L) (10) 



= f 9 (Z/L) (11) 

T,. U A 

A 

V 

where the scaling velocity is U., { = (-u’w* ) 2 , the scaling tem- 

_ w * r ' 

perature is T. ;s = — y— — , the Monin Obukhov stability length 
is L = - T o U...^/kg w ' T 1 and f^(Z/L) and f 2 (Z/L) are empirical- 
ly determined functions. 

. . 2 

An empirical expression relating to boundary fluxes 

arises from direct substitution of Equations (10) and (11) 
into Equation (8) as 

2 2 -2/3 

C T Z = T/ z A,i f 3 (Z/L) (12) 

where f 3 (Z/L) is a combination of f^(Z/L) and f 2 (Z/L). The 
formulation of f^CZ/L) was examined by Wyngaard , et al . (1971) 

and their results are shown in Figure 1. They obtained the 
following forms for f^CZ/L) 

f 3 ( Z/L) = 4.9 [1—7 (Z/L) ] — 2/3 o > Z/L (13) 

f 3 ( Z/L) = 4.9 [1+2 . 75 ( Z/L) ] 0 < Z/L (14) 

where negative and positive Z/L values correspond to unstable 
and stable conditions respectively. 
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An empirical expression relating e to boundary fluxes is 

Z Z 

Equation (10), where $^(^0 “ L can re Pl ace f^(Z/L), viz. 



eZK 

u ., 3 



4> m (Z/L) - Z/L 



(15) 



where ^^(Z/L) is the non-dimensional wind shear. Overwater 
results by Garrett (1972) on the non-dimensional dissipation 

_3 

rate, KeZU... , versus Z/L appear in Figure 2. 

Both Equations (12) and (15) can be rewritten so that 

2 

gradient terms appear in the non-dimensional C,p and e terms 
if the following expressions for T... and U... are used 



^ _ KZ 30 

* ■ 4> h (Z/L) 3 z 

TT KZ 3u 

u * ' 4 > m (Z/L) Tz 



(16) 

(17) 
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Fig. 2. 



The dimensionless dissipation rate versus Z/L. 
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Substituting Equation (16) into Equation (12) and Equation 

(17) into Equation (15) yield such expressions for the non- 
2 

dimensional C T and e terms respectively 



K' 



7 4/3,ae.2 

L 






f 3 ( Z/L) 



(18) 



£ _ r A ( Z ^ — 2 Z A ,Z,-3, 

■ 2-3u.il - K ‘♦M ( r ) • L L 1 

1 3z 



(19) 



F. RELATING Z/L TO RICHARDSON NUMBER 

The above expressions can be interpreted more readily by 
relating Z/L to the Richardson number, 



Ri = 



T( H )2 



as 



Z/L = f 4 ( Ri ) 



( 20 ) 



where fg(Ri) is an empirical function which has been well 
documented for the constant flux layer overland (Businger 
et al. , 1971) and 9 v is the vertical potential temperature, 



9 v = 0(1 + . 61q) . 



The following relationships between Z/L and Ri were 
determined by Dyer and Hicks (1970) for the unstable case 
and by Webb (1970) for the stable case: 

Z/L < 0 
Z/L > 0 



Z/L = Ri 
Ri 



( 21 ) 



Z/L = 



1-aRi 
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Here a is an empirically derived constant equal to 5. These 
relationships have been recently verified for the overwater 
regime by Paulson, et al . (1972). 

Combining Equation (18) and Equation (20) yields an ex- 

2 . 

pression for C ^ in terms of gradient measurements, 




147179072 
Z ( Tz } 



f 6 (Ri) 



( 22 ) 



The function f c (Ri) which has been determined for the overland 
6 

constant flux layer is shown in Figure 3 (Wyngaard, et al . , 
1971) . 




Fig. 3. The dimensionless temperature-structure parameter 
parameter versus the Richardson number. 
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Similarly, combining Equation (21) and Equation (19) 



yields the following gradient relationships for e 



Z 2 ( |U ) 3 
o Z 



£ 



K 2 [(l-15Ri)^ -Ri (l-15Ri) 3/4 ] Ri < 0 



(23) 




k 2 [{ 1 * 4.7 




Ri > 0 



(24) 



The curve defined by these relationships will be presented 
later (Figure 22). 

G. BULK AERODYNAMIC FORMULA 

An alternate approach for indirect estimates of the tur- 
bulent fluxes in overwater observational experiments is the 
use of bulk aerodynamic formulae. A bulk aerodynamic formula 
is one which relates the boundary fluxes to the wind speed at 
a level and the temperature and humidity differences between 
that level and the surface. Such formulae require several 
assumptions regarding the stability conditions of the boundary 
layer and the turbulent processes within it. 
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2 

Bulk aerodynamic formulae for momentum (U.,. ), sensible 
heat ( wT 7 ) and humidity (wq ' ) fluxes have the following 
forms , 




w'0 = C„ U(T -T ) = C u UAT 
Has H 

= C E U(Q a -Q s ) = C E OAQ 



(25) 

(26) 
(27) 



where the symbols have been defined on page 11 . 

Friehe (1976) summarized considerable amounts of recent 
overwater data to obtain overwater estimates of C^, C^, and 
C E « He observed and C E to be dependent on the values of 
OAT as well as height z. Friehe used the bulk aerodynamic 
formulae to formulate an alternate expression for the estima- 
tion of Z/L in the atmospheric boundary layer above the sea. 
The formulation consisted of replacing the transfer terms 
(T.,. and U... ) in the expression for L by bulk parameters 

(U A 2 = CpU, 2 and T A ) to obtain 
•* DIO 



-.231 [-^5- + 

0 


0.212 


zAQ 

- 2 J 

U 


-20 < Oat < 


25(mKs _1 ) 










(28) 


_.371[^X + 
U 


0.132 


zAQ 
-2 J 


Oat > 


2 5 (mKs -1 ) 



Equation (28) leads to a new stability parameter X, defined 
aS X = -^(1+0.212 |S) -20 < UAT < 25(mKs" ) 



X 1 



U 

^ 7 ( 1 + 0.132 j|) OAT > 2 5 (mKs - 1 ) 



(29) 
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Note that the signs of X and X* are the same as 
Z/L. 



Friehe defined X and X* for the 10 meter level. However, 
these can be generalized to any level by the approximation 

L ( ^ ) 

C n (z) C H (z) C-(z) Vz ' 0 

D _ _H _ _E _ . o . 2 nn . 

c n (io) " c„(io) " cTTToT L T ,z , J v ; 

D H E L N ( — } 

o 

This was done in this study. 

Substituting Equations (28) and (29) into Equation (12) 
yields for -20 < UAT < 25(mKs _1 ) 

22/3 -3 - 2/3 

C T Z Z z/l3 3.12 x 10 J ( 1+1. 62X) Z/J unstable 

(AT) 2 3.12 x 10 -3 ( 1-0 . 635X) stable 

for UAT > 2 5 (mKs -1 ) 

— ^ = 8.03 x 10" 3 (1 + 2.60X' )" 2/3 (31) 

(AT) 
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III. EXPERIMENT 



A. PLATFORM AND LOCATION 

Observations aboard the R/V Acania during the experiment 
in September and October 1976 were made in the area from 32° 

15 ' N to 34°13 ' N and 123°18'W to 120°20'W; Figure 4. This 
area provided a good sample for open ocean condition. 

B. INSTRUMENTATION 

Ten minute averages of temperature, wind, and relative 
humidity were taken at four levels, on two separate masts 
spatially separated on the forward deck of the ship, Figure 5. 

1. Recording 

Data logging during the experiments was accomplished 
using the NPS developed MIDAS (Microprogrammable Integrated 
Data Acquisition System) . This micro-processor based data 
acquisition system utilizes an Intel 8008 central processor 
to control the sampling, averaging, and recording of mean 
meteorological data. All software is written in PL/M to 
facilitate the writing of self-documenting program. 

2 . Mean Measurement Systems 

Mean wind measurements were made with a Thornthwaite 
Associates cup anemometer wind profile register system, model 
number 104, Figure 6. In operation, the shaft of a three cup 
anemometer unit serves as the shutter between a light source 
and a photocell for each revolution. These three plastic 
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Fig. 4. Experimental area of R/V Acania. 



legend 

Cl Cup Anemometer 
T Quart* Thermometer 
3 Humidiometer 
u* Hot Wire 
T‘ Platinum Wire 



Level 4 



Lds*il 3 



Level 2 



Level 1 




Fig. 5. Mounting arrangements aboard the Acania. 
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Fig. 6. Thornthwaite Associates Cup anemometer 
wind profile register system. 
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cone cups are reinforced with aluminum frames and are mounted 
at 120° intervals on stainless steel tubes attached to the 
main shaft. A characteristic of these cup anemometers is low 
starting speeds with a small amount of internal friction which 
aids in checking inertial overshoot. 

Mean temperature measurements were made with tempera- 
ture sensitive quartz crystal probes, (Hewlett Packard model 
HP- 2 850) at each level. (Figure 7) RF signals from the 
crystal probes and from a reference oscillator were mixed in 
the HP-2801A readout unit to produce a beat frequency whose 
signature can be analyzed to within 0.001 degrees centigrade 
per hertz. Each sensor simultaneously received pre-experiment 
calibration against a platinum resistance wire thermometer 
in a temperature controlled circulating water bath over the 
expected temperature range. The accuracy in achieving a 0.005 
degree centigrade correction factor was a constant for each 
probe. A 3.7 meter flexible coaxial cable is permanently 
attached to the sensor head and the mast mounted probes are 
housed in an aspirated shelter as seen in Figure 8. Tempera- 
ture values were automatically recorded on a printed tape. 

Mean relative humidity was obtained using the Dunmore- 
lithium chloride sensor and Hygrodynamics Digital I hygrometer 
indicator. (Figure 9) The sensor was also placed in the 
aspirated shelter. The basis of operations of this system is 
resistance charge within an electrolytic solution generating 
a reference voltage variance proportional to the relative 
humidity change. Automatic temperature compensation in the 
instruments meet the following specifications 
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Fig. 7. Hewlett Packard model HP- 2 850 temperature 
sensitive quartz crystal probes. 




Fig. 8. Aspirated shelter. 



Fig. 9. Dunmore type 

lithium chloride 
humidity sensor. 
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+_ 3% relative humidity below 90% 

+ 4% relative humidity below 90% 

Calibration of sensors was accomplished by a comparative 
method using a saturated saline solution in a closed con- 
tainer. Observation and recording was accomplished in like 
manner as the temperature. 

3 . Fluctuation Measurements 

The temperature fluctuations were measured using a 
bridge developed at GTE Sylvania, the GTE Sylvania Model 140. 
It was necessary to modify the system slightly for use 
aboard the R/V Acania. 

This system is a small temperature difference sensing 
device originally designed to be carried aloft by a balloon 
and equipped to process and transmit, by pulse-rate modula- 
tion on a radio frequency carrier, this information to a 
ground-base receiving station for demodulation and recording. 

The baseband portion of this system is a balanced 
Wheatstone bridge excited by a 3KHZ signal with a synchronous 
detector on the output. Segments of a very small diameter 
platinum wire serve as temperature sensors in opposite arms 
of the bridge. In the single wire mode, one wire is replaced 
by fixed resistance. An output from the bridge which is 
proportional to the temperature difference between the two 
probes or a temperature change of a single wire in the single 
wire mode is produced. 

— 6 

The sensor wire is 0.5 cm long and 2.5 x 10 m in 
diameter. This extremely small mass allows a response to 
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temperature variations of up to 1 KHZ while electronic ampli- 
fication allows temperature difference as small as 0.004 
degrees C to be observed. 

Wind fluctuations were recorded on a Honeywell model 
5200 14 channel tape recorder. Real time readout on an eight 
channel chart recorder. Brush model 240, was used to check 
the quality of the signals coming from the sensors. 

4 . Fluctuation Statistics 

Spectral analyses to obtain values were performed on 
fluctuating velocity data obtained with Hot Film sensors 
oriented in the vertical, normal to the mean wind. These 
data represent a point measurement and yield temporal or 
time descriptions of the fluctuations. Ten minute segments 
of data originally recorded on magnetic tape were recorded 
into an EMR-Schlumber model 1510 digital spectrum analyzer. 
Procedures for converting spectral values, obtained with the 
analog spectral analyzer, to engineering units and for obtain- 
ing turbulence parameters from the spectra are described in 
the following paragraphs: 

a. Scaling Spectral Plots 

A necessary procedure was scaling the spectral 
plots to relate RMS input voltages to power spectral densi- 
ties (PSD); variance per unit frequency. To obtain power 
spectral density levels, corresponding to RMS voltage inputs, 
calibrated scale charts had to be constructed. 

For purposes of the x-y plot format of the 
analyzer output, the RMS voltages were converted to y = log. , 
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(voltage) units and a graduated scale was constructed so 
that the logarithm of volts RMS could be interpolated from 
spectral plots. The value of the vertical scale (y) was 
adjusted for each spectrum as a function of both input gain 
and spectral gain. These values were then converted to PSD 
levels for use in calculating e values. 

The relation between Volts RMS to PSD units used 



was 



S(f) PSD units 



(cal. level VRMS) 2 /n „y N 2 
1.5 Bandwidth x '‘ 1U ' 



(32) 



where 



Bandwidth 



Freq. Range 
No. Channels 



( 



256 Hz 
Hz 



1 Hz) 



and cal. level VRMS 

( 



voltage at y = 0 

1 VRMS for 3.16 V input setting). 



Amplitude scaling calibrations were accomplished 
using externally generated "white noise" signals of 1 Volt 
RMS with a frequency range from 0.1 Hz to 1000 Hz (giving a 
PSD of 10 -2 V 2 /Hz). Setting a 3.16 V (lOdB) input on the EMR 
1510 Digital Spectrum analyzer insures that 1 Volt corre- 
sponds to y = 0. An example of such a calibration plot for 
a frequency range of 1-200 Hz (the frequency range of inter- 
est) is shown in Figure 10. 

b. Computation of Turbulence Parameters 
from Scaled Spectra 

Values of the turbulence parameter £ were obtained 
from the velocity variance spectra on the basis of the formula 
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Fig. 10. Spectrum calibration plot. 



for the inertial subrange in wave number space S(k), Equa- 
tion (6), which predicts a -5/3 slope for the spectra when 
plotted in a log-log format. Figure 11 is a typical spectra 
considered in the analyses. Velocity variance spectra con- 
sistently exhibited slopes of -5/3. 

Assuming -5/3 slopes for the variance spectra, 
the intercept of the best fit -5/3 slope with the 1 Hz fre- 
quency line was the spectral density denoted (PSD) value used 
in computing the parameter of interest, e . The measured PSD 
value was converted to engineering spectral density units by 
the relation 

S u (f) = C u 2 • PSD 

= ( cm/s/Volt ) 2 • Volt 2 /Hz =(cm/s) 2 /Hz (33) 

where C u was the calibration factor for the constant tempera- 
ture, anemometer. Determination of the value for has 
been described by Lund (1975). 

Since the velocity fluctuations were measured at 
a fixed point in the flow, the resultant spectral values are 
defined at "temporal" frequencies denoted as S (f) in Equa- 
tion (33). 

5 . Analysis of Paired Temperature Sensor 

2 

C^, values were estimated by using paired platinum 
wire sensors for variance analysis. This analysis was based 
on the expression for the temperature structure function, 
Equation (5). The analysis is one in which temperature differ- 
ences between the paired resistance wires separated by a 
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Fig. 11. Velocity spectrum. 



distance r is subjected to variance analyses as indicated by 

the term (T'(x) - T'(x+r)) 2 . 

Voltages corresponding to the temperature differences, 

AT, as measured by the paired sensors were recorded on analog 

magnetic tape. The variance or root mean square of the tem- 

2 

perature differences required to compute was obtained by 

analyzing the analog voltages with a mean square voltmeter 
CThermo Systems Incorporated, model 1060) using a 30 second 
time constant for the averaging. The output of the mean-square 
voltmeter (as well as a record of the original signal) were 
recorded on standard strip charts. 

Ten minute averaged values of the RMS voltages were 
abstracted from the strip chart and these values were utilized 
in a program for the IBM 360 computer. This program con- 
verted the mean square voltages values to temperature differ- 

2 

ence variance values, from which were computed. 

6 . Fluctuating Sensor Calibration 

Platinum wire resistance is related to temperature by 

R = R (1 + cc(T-T )) ( 3 U) 

o o 

where R = resistance, T = calibration temperature, T q = 0°C, 

R q = resistance at T q , and a = .0036/°C (temperature coeffi- 
cient of resistance for platinum). Each wire's resistance 
was measured in the laboratory at a temperature T, to deter- 
mine R . 
o 
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Temperature fluctuations can be related to wire re- 



sistance fluctuations by differentiating Equation (34) 

dR = ctR dT 
o 

or 

R' = aR Q T ' (35) 

Recorded voltage fluctuations are then related to re- 
sistance changes by the transfer function, H, in the expres- 
sion 



V ? = HR' 



(36) 



where H was determined for the bridges by varying the resis- 
tance of a TSI 1056 Variable Decade Probe Resistance and 
recording the corresponding voltage changes. Voltage vs 
resistance is plotted and the slope of this plot determines 
the transfer function, H. For this experiment H = 10 volts/ 
ohms . 

Finally, to relate recorded voltage fluctuations to 
temperature fluctuations, Equation (36) is substituted into 
Equation (35) giving 

T* = ( 1/H R q ) V' (37) 

or 

T' = C H V' (38) 

where is the temperature calibration factor in °C/volt. 



40 



Velocity sensors calibrations were accomplished in a 
wind tunnel using the Hot Wire Calibrator TSI 1125. In this 
process, a known wind velocity is blown across the wire and 
the corresponding voltage output is recorded. The sensor wind 
speed is given by 

V 2 = aU 1 ' 5 + b (39) 

where V is the hot-film voltage output and 0 is the mean wind 
speed for any given level. The constants, a and b, are the 
calibration curve slope and intercept respectively. The con- 
stants, a and b, were derived in the lab. 

To convert the voltage PSD levels to velocity units 
requires a calibration factor given by 

U T = CV* (40) 

where C is the calibration factor in m/s/volt, V* is the 
voltage fluctuation, and U T the velocity fluctuation. 
Differentiating Equation (39) yields 

U* = (4vtf*/a) V' (41) 

which when substituted into Equation (40) yields 

C = 4VU ; Va (42) 

The slope of the plotted curve yields the value for a in the 
above equation. 
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C. ANALYSIS PROCEDURES 



The data were edited for gross errors or inconsistencies 
due to known instrumentation malfunctions or limitations of 
instruments. An example of limitation on instruments is the 
relative wind direction being greater than 30° off the bow of 
the ship. Earlier experiments have shown that when these con- 
ditions exist, the winds are greatly influenced by the ship 
and therefore are of little value to the experiment. 

The criteria for evaluating or discarding data from in- 
dividual levels or for entire interval depended on performance 
check results obtained during measurements, obvious inconsis- 
tencies, and the availability of fluctuation data. Since mean 
wind 0 , mean temperature T, and mean humidity a, are expected 
to vary logarithmically with height, these parameters were 
plotted for each time period using the IBM 360 and NFS FORTRAN 
program PLOTP. A best fit straight line was drawn to the data 
points. (See Figure 12.) This procedure was very subjective, 
but was intended to further refine the data by detecting more 
obvious discrepancies and determining if the data did indeed 
vary logarithmically with height. Gradient values of the 
parameters used in actual computations were evaluated using a 
least square fit. By means of the subjective procedure, it 
was determined that the humidity sensor at the second level 
was reading too low during most of the experiment. (See Figure 
11a.) This error was corrected for the least square fit by 
assigning the second level the same value as level three. 
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Fig. 12. Typical profiles for a) specific humidity - 25 Sept 
75 (0105); b) temperature - 25 Sept 75 (0250); and 
c) winds - 25 Sept 75 (0551). 
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IV. RESULTS 



A. C T 2 RESULTS 

1. Gradient Description 

2 

A comparison of 36 8 overwater Crp and profile results 
with the overland empirical expression, Equation (22), appears 
in Figure 13. The distribution predicted by Equation (22) is 
the solid line. The circles are mean dimensionless tempera- 
ture structure function parameter (DTSFP) results over .50 
Richardson number (Ri) intervals. The Standard deviation 
corresponding to each mean value appears as vertical lines 
on which the number of values defining the mean appears. Sig- 
nificantly, these results represent more than a factor of 3 
(368 versus 119) of those presented previously. (Davidson, 
et al., 1976). However, the mean DTSFP distribution changed 
very little between the two sets. Larger DTSFP values for 
the overwater regime for near neutral and stable conditions 
(Ri > -.5) can be qualitatively related to the wave influence. 
Decreased DTSFP values for unstable conditions (Ri < -.5) are, 
perhaps, associated with aerodynamic smooth properties for 
the sea surface during light wind conditions. 

A curve describing the observed data was obtained 
using a best fit polynomial analysis. The results of this 
analysis were 
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= 2 . 37( . 74+Ri) -5 • 04 Ri > 0 (43) 

C 2 /Z 4/3 (|i) 2 = 0 - 56 ( . 24 + Ri)' 2 * 13 -1.25 < Ri < 0 

1 d Z 

(44) 

=0.2 Ri < -1.25 (45) 




Fig. 13. Overwater results for DTSFP versus Richardson 
number . 

The influence of the waves can be further examined on 

2 

the basis of observed height variations of Crp . Similarity 

2 . ... 

theory predicts that C^, has a height variation which depends 

on stability which was discussed by Wyngaard et al. in con- 

2 . 

junction with an examination of Equation (12). Since T... is 

independent of height in the surface layer, Equation (12) 

yields the following general predictions for the variation of 
2 

C,p with height based on stability. 

2 

a. For near neutral conditions, (Z/L = 0), de- 

-2/3 

creases with height as Z Equation (12) becomes 

2 2 -2/3 

C T ^ = T ... 1 Z 6 (46) 
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b. For unstable conditions (Z/L <_ 0, Equation (13)), 

2 -2/3 

Crp decreases more rapidly with height than Z 

2 

If -7Z/L >> 1, C^, would decrease with height as 
-4/3 

Z , since Equation (12) becomes 

C T 2 ~ 4/3 j~- (wT^) 3 Z“ 4/3 (47) 

with T 2,/3 /kg(wT' ) 4/, ~-' independent of height. 

2 

c. For stable conditions (Z/L _> 0),C^, will decrease 

-2/3 

with height less rapidly than Z since 

C T 2 = T.,. 4. 9Z -27 3 [1+2.75 |] (48) 

and the term in brackets increases with increas- 
ing height. 

The wave influence can also be seen in height versus 
2 2 2 

C T or C N (computed from C T values) results in Figures 14 

and 15 respectively. In Figure 14, the feature which suggests 

-2/3 

the wave influence is the departure from the Z height dis- 

tribution as the surface is approached. In the upper three 
2 -2/3 

levels Cfj, appears to fit a Z distribution, but as the 

surface is approached the distribution appears to be better 

-4/3 

defined by a Z height distribution. In Figure 15, a 

similar change is observed. The upper levels appear to fit a 
_ 2 / 3 2 

Z distribution for C N but as the surface is approached 

. o 

especially at the one meter level, the values appear to 

-4/3 

fit a Z height distribution. 

2 

Further evidence of wave influence on appears in 

joint probability-conditional mean results of the normalized 
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Fig. 14. 



. 2 

Normalized C T results versus height (Z), 114 
periods . 




Fig. 15. Normalized (computed from C T 2 ) results 

versus height. 1 
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temperature variance (a^/T... ) over waves in Figure 16. The 

results and the measurements yielding them have been de- 

2 2 

scribed by Davidson (1974). The ratio (a,p /T.. ; ) is direct- 
ly related to the DTSFP, defined by Wyngaard, et al. (1971). 
In Figure 16, the vertical axis variable (C/U...) defines wave 
influence on cr^/T... , C is the phase speed of prevailing swell, 
L is the Monin-Obukhov stability length. The results were 
obtained under stable conditions. The diagonal orientation 
of the contours of the number joint occurrences provides 
qualitative evidence that the waves' influence (C/U... ) on 
o t /T,., hence the DTSFP, is as significant as the stability 
influence (Z/L). 




Fig. 16. Joint probability-conditional mean results of 
a m/T... for C/U... and Z/L dependence. 
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2 . 



Bulk Aerodynamic 



Results obtained from analyses with bulk aerodynamic 
variables, Equation (31), are compared to the Wyngaard et al. 
prediction curve in Figure 17. The latter is the solid line. 
The circles are mean dimensionless structure function para- 
meter results over 1.00 X intervals. Standard deviations 
corresponding to mean values appear as vertical lines and the 
number above these vertical lines indicates the number of 
values defining each mean. Larger values of DTSFP are again 
attributed to the wave influence as explained earlier. 

The following was evaluated for the bulk aerodynamic 
expressions using a best fit polynomial analysis. 

= . 0 3 7 ( . 642 + X)" 1, 72 Stable (49) 

c t 2 z 2/3 /(at) 2 

= . 042 ( . 3 2 3-X) “ * 5 7 Unstable (50) 

Three points from an analysis by Friehe (1976) appear 

in Figure 17. A significant aspect of Friehe' s analysis of 

2 

available data was the exclusion of C^, data from temperature 
fluctuation records which contained large numbers of spikes 
exhibiting ramplike signatures, as shown in Figure 18. This 
led to a significant reduction in the amount of data used in 
his comparison. Based on observations by Hicks (1972) and 
hypothesis of Schmitt et al. (1976), Friehe attributed the 
occurrence of the spikes to salt nuclei accumulations on the 
wires with subsequential erroneously sensed temperature 
changes by the system as water vapor is absorbed or is 
evaporated on the salt nuclei. 
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results for dimensionless temperature-structure parameter 




Fig. 18. Time series showing "cold spikes" in the 

temperature field under stable conditions. 
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The occurrence and nature of the spikes were also 



examined with data collected over Lake Michigan by Davidson 
and Safley (1972). Figure 19 includes ramp signatures in 
temperature traces observed over Lake Michigan where salt 
accumulation problem would not be a factor. Therefore, the 
ramp signatures is possibly due to factors associated with 
the nature of the shear flow itself. On this basis there 
did not appear to be sufficient evidence to discard data 
with spikes with ramp signatures. It is not believed that 
this data being included in this study caused the signifi- 
cant difference between the predicted Wyngaard curve and the 
observed data curve. 

A contradiction occurs between results in Figures 13 

9 0 

and 17. In Figure 13, given the same Z and values for 

both regimes, the plotted curves predict a higher value of 

2 ... 

Crp for the overwater regime for Ri > -.5 and a higher value 

2 

of C^, in the overland regime for Ri < - . 5 . On the other 

hand in Figure 17 , given the same Z and AT values for both 

2 . 

regimes, the plotted curves predict a higher value of C m m 

2 

the overwater regime for X < 1.5 and a higher value of C_ 
in the overland regime for X > 1.5. Since Ri and X are cor- 
related they should predict basically the same result. How- 
ever, as is seen in Figures 13 and 17, the only ranges where 
they correlate is the near neutral range . 

In order to explain the reason for this contradic- 
tion, values of Ri were plotted versus values of X, Figure 20. 
The circles show average values of Ri for X intervals of • 
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Fig. 19. Temperature trace and velocity fluctuation vectors 
and wave heights associated with microthermals 
observed over Lake Michigan. 
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Fig. 20. X versus Richardson 
number. 
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Since Z/L is correlated with X, a similar relationship 
should be seen between X and Ri as exists between Z/L and Ri . 
The depicted curve is similar to the latter relationship. 

The conclusion is that the expected correlation existed be- 
tween these two stability parameters. 

A second possible reason for the contradiction is the 

relationship between Z^^/AT^and l/Z^ 7 ^ ( 30/3z) ^ which norraa- 

. 2 . 

lize Crp in the DTSFP. Mean values of the former versus the 
latter appear in Figure 21. There should be a positive cor- 
relation between the two and this occurs only for large values 
of l/Z^ ^ 09 / 3z ) ^ . For small values of the latter there is, 
however, a negative correlation. This result indicates that 
the two parameters are correctly correlated when small tem- 
perature gradients exist, but are not when the gradients be- 
come large. The contradiction is likely to be due to the 
■ problem of accurately defining the bulk temperature difference, 
AT. 

B. e RESULTS 

1 . Gradient Description 

Comparison of the observed £ results and the curve 
based on Equations (23) and (24) is shown in Figure 22 where 
the predicted curve is the solid line. Averages over Ri in- 
tervals of .25 of observed data appear as small circles. 

The error bars are standard deviations from the mean within 
each interval. The number of observations defining each 
mean value is at the top of the error bars. The significant 
difference between the predicted curve and the observed data 



55 




CM 

✓ s 

O |N) 

co 



CN 



N 



DC 



J 5 

fO^ 


cn 


jn 


s 


N 


c 




> 



CM 



ESI 



CM 



W3 

•H 

U* 



4? 

CN 

N 



CN 



56 




57 



The dimensionless dissipation rate versus Richardson number. 



for Ri > .20 is expected on the basis of Equation (23). For 

Ri > .20 the equation is no longer valid since terms in the 

expression become infinite. 

An examination of Equations (13), (14), and (15) 

3 

with the knowledge that U... and L are independent of height 
in the surface layer yields the following general prediction 
for the variation of e with height: 

a. For neutral conditions, (Z/L =0), e will decrease 
with height as Z - "* - since Equation (15) becomes 




b. 



c . 



For stable conditions, (Z/L > 0), s will decrease 
less rapidly with height than Z ^ since Equation 
(15) becomes 



U,.' 



£ = 



K 



Z -1 ( 1+ 3 . 7 Z/L) 



(52) 



For unstable conditions, (Z/L < 0), e will de- 
crease less rapidly with height than Z ~ since 



Equation (15) becomes 

3 



£ = 



U ... ' 

K~ 



Z -1 ( 1-15 Z/L)" 5 * -Z/L 



(53) 



The quantity in brackets increases with increas- 
ing height. This increase is less than for 
stable conditions. So the decrease of £ with 
height for unstable conditions is less than that 
for stable conditions. 
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Results in Figures 23 and 24 indicate that under near neutral 
conditions e fits a Z ^ height distribution quite well. 

There appears to be no change in the height distribution as 

. 2 
the surface is approached as there was with . 

2 . Bulk Aerodynamic 

Results obtained from analyses with bulk aerodynamic 
variables appear in Figure 25. The predicted curve is the 
solid line, while the circles are mean non-dimensional e re- 
sults over X intervals of 1.0. Again, standard deviations 
corresponding to mean values are vertical lines and the 
number above these lines indicate the number of values defin- 
ing each mean. The correlation between the non-dimensional 
s values versus X, as described for Z/L by Busch and Panofsky 
(1969), is quite good. 

C. C T 2 RESULTS AT REGIONS ABOVE SURFACE LAYER 

A Cessna was used in this experiment to extend the data 

base beyond the height of the instrumented towers. Height 

2 

distribution of C^ values as expressed in Equations (46), 
(47), and (45), have been shown to be valid to height above 
4000 meters (Tsvang, 1969 ). Figure 26 shows that a Z 
distribution exists, indicating stable or near neutral condi- 
tions, from lowest level to about 100 meters. Above 100 

2 

meters, there is a rapid increase of C,p values. 

Frisch and Ochs (1975) suggested that if the height of 
the inversion base, Zi, is another relevant scaling parameter, 
then proceeding in a manner similar to that of Wyngaard et al , 
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Figo 23. e versus log Z. 
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Fig. 24. e versus log Z 



CO I— GH 



i 




x 



LO 

CM 

3 

• r-j 

u« 



62 



The dimensionless dissipation rate versus 







T(*C) 



Fig . 26 . 



Crp height distribution and temperatur 



profile . 
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Equation (47) could be modified with an additional non-dimen- 
sional function, G(Z/Z^). Equation (47) then becomes 

C 2 = [ — ~] (T/g) 2 / 3 (w r 9 T ) 4 / 3 Z _4/3 ) (54) 

1 3k Z/J Z i 

where G(Z/Z^) = 1 for Z/Z^ << 1 in order to have Equation 

(54) reduce to Equation (47) in the surface layer. 

. . . . _ u / 3 

Frisch and Ochs found a significant departure from a Z 

law for values of Z/Z^ > 0.1. This indicates that the atmos- 
phere is no longer in the "free convection" regime, but is 
affected by the lid on the convective layer at Z^. This 
effect would restrict the integral scale of the turbulent 
temperature fluctuations with height compared to the free con- 
vection value. In addition, because of the layer scale 
circulation and mixing, the temperature variance should be 

greater than the free convection variance. These effects 

2 

will increase the value of C,p relative to the free convec- 

2 

tion prediction. Figure 26 shows the rapid increase of 
values as the inversion is approached. 
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V. CONCLUSIONS 



Application of formulated overland expressions in the 

overwater regime showed a significant difference between ob- 
2 

served C^, overwater data and overland predictions . This 
difference is attributed to the wave influence for near 
neutral and stable conditions and to the aerodynamic smooth 
properties of the sea surface during light wind periods for 
moderate unstable conditions. The formulation of these over- 
water curves will greatly enhance the Navy's operational 
ability to describe the marine boundary layer and its in- 
fluence on the optical wave propagation. 

In comparing the DTSFP overland and overwater when plotted 
against Ri and then versus X, significant differences were 
noted. Conclusive explanations for this difference cannot 
be given in this study. Therefore, recommend that serious 
consideration be given to an in-depth examination of this 
significant feature. 

Correlation of non-dimensional e values with the stability 
parameter, Ri , showed agreement between results and predic- 
tions. For Ri > .20, the theory is not valid and therefore 
there was very poor correlation. Any attempt to describe 
the normalized e for Ri > .20 in the future will require an 

expression which will remain valid for all values of Ri . 

2 • . 

Crp results appear to be described by the same high dis- 
tributions up to levels well above the surface layer, 
approximately 60% of the distance to the marine inversion. 
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